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We demonstrate that parity nonconserving interaction due to the nuclear weak charge Qw leads 
to nonlinear magnetoelectric effect in centrosymmetric paramagnetic crystals. It is shown that the 
effect exists only in crystals with special symmetry axis k. Kinematically, the correlation (correction 
to energy) has the form Hpnc oc QwE ■ [B x k] (B ■ k) , where B and E are external magnetic and elec- 
tric fields. This gives rise to the magnetic induction Mpnc oc Qw {k(B ■ [k x E]) + [k x E](B ■ k)}. 
To be specific, we consider rare-earth trifluorides and, in particular, dysprosium trifluoride which 
looks the most suitable for experiment. We estimate the optimal temperature for the experiment 
to be of a few kelvin. For the magnetic field B = 1 T and the electric field E = 10 kV/cm, the 
expected magnetic induction is 47tMpnc ~ 0.5 • 10 -11 G, six orders of magnitude larger than the 
best sensitivity currently under discussion. Dysprosium has several stable isotopes, and so, compar- 
ison of the effects for different isotopes provides possibility for precise measurement of the Weinberg 
angle. 
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I. INTRODUCTION 

Studies of atomic parity nonconservation (PNC) pro- 
vide important tests of the standard model of electroweak 
interactions and impose stringent constraints on the new 
physics [| . The effects of PNC in atoms due to the nu- 
clear weak charge Qwhave been successfully measured 
for bismuth 0, lead thallium Q, and cesium 
The authors of 0, have also performed measurements 
of the nuclear spin-dependent PNC effect in Cs caused 
by the nuclear anapole moment (e.g., see text Q), which 
remains the only measurement of this kind to date. 

In the present work we concentrate on effects caused by 
the nuclear weak charge in solids. We predict a new kind 
of (nonlinear) magnetoelectric effect in centrosymmetric 
crystals due to the violation of parity (P) at fundamental 
level, proportional to Qw- Although we show that the 
experimental observation of the proposed effect is pos- 
sible at the current level, the uncertainty of theoretical 
calculation would not allow for an accurate interpretation 
in terms of the nuclear weak charge Qw itself. Instead, 
the potential importance for fundamental physics lies in 
the isotope dependence of the effect we suggest. We ar- 
gue that comparing the value of the effect in samples with 
different isotope content, one can, in principle, determine 
the Weinberg angle with high precision. In this respect 
the idea is similar to the suggestion presented in Ref. Q 
of an atomic PNC measurements in rare-earth Dy, and 
the later suggestion of Ref. 9] for a measurement with 
atomic Yb. 

The idea to use solids for studies of fundamental sym- 
metry violations dates back to 1968, when Shapiro sug- 
gested that the existence of the permanent electric dipole 
moment (EDM) of electron should lead to linear magne- 
toelectric effect in solids 10]. (The existence of the per- 
manent EDM of a quantum particle violates both parity 
(P) and time-reversal (T) symmetries, see, e.g., [Tlj.') 



However, the results of the early experiment 0] (1978) 
performed in Ni-Zn ferrite were not impressive due to 
experimental limitations. In the recent paper |l3j . Lam- 
oreaux suggested that application of novel experimental 
techniques in garnet materials should lead to substan- 
tial increase in sensitivity compared to the previous at- 
tempt; this suggestion renewed the interest in the original 
idea. Calculations 0] have shown that improvement of 
the statistical sensitivity to the electron EDM of several 
orders of magnitude was possible compared to the cur- 
rent experimental limit on this value |l5|. At the same 
time, the experimental search was initiated — the first re- 
sults are already available in the literature, 0] (see also 
Ref. [13). The other suggestions for fundamental sym- 
metry violation search in solids include proposed search 
for manifestations of the nuclear anapole moment in gar- 
nets 0] (violation of parity), and the search for the 
nuclear Schiff moment in the ferroelectric lead titanate 
compound 0] (violation of T and P symmetries). 

One of the most important issues for any solid state ex- 
periment is certainly the existence of various systematic 
effects. It is getting more and more clear that the right 
choice of compound is crucial. In particular, any kind 
of spontaneous magnetic ordering in the crystal (ferro-, 
antifcrro-, or fcrrimagentism) should be avoided, thus, 
we consider paramagnetic material — rare-earth trifluo- 
ride RF3. The PNC weak interaction of electrons with 
the nuclei of rare-earths gives contribution to the total 
energy of the crystal 

% c «Q ff (E.[Bxk])(B-k) , (1) 

in the applied electric field E and magnetic field B; vec- 
tor k is directed along the crystal symmetry axis. Cor- 
relation similar to Eq. Q> was considered previously by 
Bouchiat and Bouchiat 20] in regard to their proposal to 
measure nuclear anapole moment using Cs atoms trapped 
in solid 4 He. 
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The effective interaction (QJ induces the following PNC 
macroscopic magnetization which can be measured in a 
magnetometry experiment: 

M = oc Qw{HB ■ [k x E]) 

+ [kxE](B-k)} . (2) 

Naively, one can expect the PNC magnetization of the 
form M oc [B x E], however, this kind of magnetization 
does not lead to any energy shift in the applied mag- 
netic field, and hence, does not appear. Eq. presents 
the simplest energy correlation allowed kinematically. To 
generate this correlation we need a crystal that allows 
nontrivial second rank material tensors, so the crystal 
symmetry must be lower than cubic. In addition, we 
need a necessarily centrosymmetric crystal in order to 
avoid the imitation of the PNC effect by spontaneous 
parity breaking in the lattice. 

Rare-earth trifluoride compounds have hexagonal 
structure described by C^.i space group and, apparently, 
satisfy the crystallographic criteria we set. These com- 
pounds are widely used as laser materials, and as such, 
their electronic structure has been studied extensively 
and is rather well-understood. As it will be seen later, 
DyF 3 is, probably, especially well-suited for the exper- 
iment because of the large value of the PNC effect in 
this compound. Dysprosium has also particularly large 
number of stable isotopes which can be used for the de- 
termination of the Weinberg angle as mentioned above. 

The structure of this paper is as follows. Section II 
deals with atomic estimates involved in the problem. 
Starting from the electron-nucleus weak interaction, we 
consider an ion in the lattice environment and derive the 
effective single-ion PNC Hamiltonian. In Section III we 
average this Hamiltonian over the crystal structure of the 
compound, which gives the kinematical structure of the 
effect and its dependence on the external fields. Section 
IV presents our results and discussion of uncertainties 
involved. 



II. SINGLE ION EFFECTIVE HAMILTONIAN 

The weak charge-induced PNC interaction between 
electron and nucleus is of the form (see, e.g., 0) 

Hw = -^=5(r) 75 Qf. (3) 

Here G is the Fermi coupling constant, <5(r) is the Dirac 
delta function for electrons representing the nuclear den- 
sity distribution, 75 is the Dirac 7-matrix. The nuclear 
weak charge is Qw = — N + Z(l — 4 sin 2 9w) where 6w 
is the Weinberg angle, N is the number of neutrons, and 
Z is the number of protons. 

Let us consider an ion embedded in crystal lattice, and 
let us impose an external electric field E which shifts 
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FIG. 1: A typical atomic perturbation theory diagram re- 
sponsible for the PNC Hamiltonian @. The wavy line on the 
diagram represents the residual Coulomb interaction between 
electrons, the weak interaction Q is shown by the cross on 
the electron line, and the dashed line denotes interaction of 
electron with the external electric field E. 

the ion from its equilibrium position with displacement 
X. The ion has the orbital angular momentum L and 
the spin S, in this section we assume that L and S are 
decoupled from each other. Certainly, in reality L and 
S are coupled due to the combined action of the spin- 
orbit interaction and the non-central crystal field of the 
lattice at the ion site. We will consider these effects in the 
next section, but for now we switch these interactions off. 
The interaction @ is a T-even pseudoscalar, therefore, 
the effective Hamiltonian which arises after calculation of 
the electronic matrix element of (01 must be of the form 

UPNC oc E • [L x S] . (4) 

This is the only structure allowed kinematically, and the 
present section is devoted to the calculation of this effec- 
tive Hamiltonian. 

Rare-earth trifluorides RF3 are ionic crystals contain- 
ing triply ionized rare-earth ions R 3+ . The general elec- 
tronic configuration of these ions is ls 2 ...5s 2 5p 6 4/ n with 
the number of 4/-electrons ranging from 1 in Ce 3+ to 13 
in Yb 3+ . So, L and S are carried by the 4/-electrons 
which constitute the only open shell, /-electrons practi- 
cally do not penetrate to the nucleus and hence, practi- 
cally do not interact with the weak charge directly. In- 
stead, the effective interaction arises in the third order of 
perturbation theory due to many-body effects, a typical 
diagram is shown in Fig. ^ The vertical wavy line on 
the diagram connecting the two electron lines represents 
the residual Coulomb interaction between electrons; the 
weak interaction J3J is shown by the cross on the electron 
line; the dashed line denotes interaction of electron with 
the external electric field E. This is but one of a number 
of diagrams responsible for the effect, in order to find the 
effective interaction (@J one has to evaluate all of the dia- 
grams. A similar atomic calculation has been performed 
previously in the paper |lSj | in regard to the interaction of 
atomic electrons with the nuclear anapole moment. The 
effective interaction with the anapole moment appears 
in the third order of perturbation theory as well, and, 
although it has different kinematical structure, the part 
involving atomic radial integrals is practically identical 
to the present case. The main difference is that the weak 
charge matrix element (Hw), Eq. shown by the cross 
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in the diagram Fig. has to be replaced with the ma- 
trix element of the anapole interaction (H a ). Therefore, 
to estimate the effective Hamiltonian (@J, it is not nec- 
essary to repeat the atomic calculations — it is sufficient 
to rescale the results of Ref. . Analytical expressions 
for both the weak charge (Hw) and the anapole moment 
(H a ) interaction matrix elements based on the semiclas- 
sical approximation for electron wave functions can be 
found in the book Jjj, and lead to the following ratio: 

(s 1/2 \H w \p 1/2 ) ^ 
(si/2\H a \pi/ 2 ) n a 

Here K a is the dimensionless constant characterizing the 
nuclear anapole moment, it is of the order n a ~ 0.4 . 
Rescaling the result of Ref. we find the the single- 
electron effective Hamiltonian for 4/-electrons: 

ffPNC~2-Hr 16 (X-[lxs])Q w Eo . (6) 

Here 1 and s are the single-electron orbital angular mo- 
mentum and spin, and Eq — 27.2 eV is the atomic energy 
unit. The ionic displacement X is measured in units of 
the Bohr radius as- 

Simple electrostatic considerations allow to relate the 
displacement X of the triply ionized rare-earth ion in the 
lattice to the strength E of the external electric field: 

P = 3enXa B = !-E ; (7) 



47r 3e naB 

Here P is the dielectric polarization, e = |e| is the elemen- 
tary charge, n f=a 1.8 • 10 22 cm~ 3 is the number density 
of rare-earth ions in the compound, and e « 14 is the 
dielectric constant. Evaluating Eq. © we get 

X[a B ] = 2.5 • ltr 8 E[V/cm] . (9) 

Finally, we have to average the single-electron effec- 
tive Hamiltonian |JBJ over the ground state of the 4/™ 
electronic configuration of the R 3+ ion with the total 
momenta L = Ii + ... + \ n and S = si + ... + s„ (see also 
the remark [21|). The averaging was performed with the 
help of the usual fractional parentage coefficients, see, 
e.g., . The result reads: 

iW ~ A • 1(T 24 E • [L x S] Q w E , (10) 

where the electric field E is measured in units of V/cm. 
Values of the coefficient A in the above equation for dif- 
ferent rare-earth ions are listed in Table|U Note, that the 
effect vanishes for Gd 3+ ion which has the ground state 
with L = 0. 

Eqs. ©, i(TU)) were obtained by rescaling the results of 
our previous calculation [Tsj . where the rare-earth ions 
were considered in the garnet environment consisting of 
eight oxygen O 2- ions. In the present case the lattice 
structure is different, and the environment is composed 
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1.7 


1.0 





-1.0 
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non-Kramers ions 




Pr 


Pm 


Eu 


Tb 


Ho 


Tm 




L 


5 


6 


3 


3 


6 


5 




S 


1 


2 


3 


3 


2 


1 




A 


2.5 


1.3 


0.8 


-0.8 


-1.3 


-2.5 





TABLE I: Values of the coefficient A in the effective parity 
nonconserving Hamiltonian llOH . For quick reference, we also 
present the values of the orbital momentum quantum number 
L and the spin S for each rare-earth ion. 



of fluorine ions F _ . For the present calculation which 
concerns the Hamiltonian IjlOfl. however, only the spher- 
ically symmetric part of the electron density from the 
neighboring ions plays role. Since the ions O 2- and F~ 
have the same electronic configuration, 2p 6 , and the in- 
teratomic separation in both compounds is roughly the 
same, we believe the rescaling procedure gives reliable 
results. At most, we can imagine an error in the coeffi- 
cient A, Eq. HI U| l. °f a factor 2 — 3, which is acceptable 
for the purposes of the present estimate. In the next Sec- 
tion, when we consider coupling to the external magnetic 
field, the anisotropic nature of the environment is very 
important. There, we treat the ionic environment more 
accurately in the framework of the Coulomb crystal field 
model. 



III. CRYSTAL STRUCTURE AND INDUCED 
MAGNETIZATION 

In the present section we take calculation of the PNC 
interaction Hamiltonian (|1U|1 one step further. The idea 
is that the PNC part is exceptionally small compared to 
the other terms in the full Hamiltonian for the ion in 
external magnetic field: 

H = A LS (L-S) + H c[ + fi B B-(L + 2S)+Hp NC . (11) 

(Here A^s is the constant for fine-structure splitting, H c f 
denotes the crystal field interaction, /ie is the Bohr mag- 
neton.) So, it is sufficient to consider the iJpNC in the 
first order of perturbation theory only, 

ASp NC = (<7|#pnc|<7) . (12) 

The ground state \g) is determined by the first three 
terms in the Hamiltonian (|11|) . 

For the case of isotropic environment, the direction 
of (L) and (S) is completely determined by the external 
magnetic field: (L), (S) cx B; the spin-orbit coupling con- 
serves J = L + S, so that (L), (S) cx (J), and (J) cx B. 
Intuitively, it is then clear that the expression (L x S) 
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in (|10fl must vanish; the result is rather easy to prove 
rigorously if one considers the ground state \g) — | J, J z ). 
So, the PNC effect (|10|) vanishes for environments with 
spherical symmetry, and the same result is expected for 
cubic lattices. The situation is different, however, for an 
ion in the crystal environment with an axis k. This is 
especially obvious in the case of zero spin-orbit coupling 
when (L x S) = (L) x (S) . Indeed, the direction of (L) is 
then determined not only by the external magnetic field 
B, but also by the lattice crystal field: in the simplest 
case we can write (L) oc B + (B • k) k. The total spin (S) 
is oriented along the magnetic field. Then, the value of 
the cross product (L x S) oc (B • k)[B x k] is, generally, 
not zero. Existence of the LS'-coupling certainly compli- 
cates the analysis, but the general picture remains the 
same. 

The crystal field Hamiltonian H c { in Eq. 1|11[1 was cho- 
sen in the form suggested in paper |23|: 



H c f = PnAnm ^ J ^ ~ 1 Y nm fa ) , (13) 

nra i 

where is the radius-vector of the 4/-electron, and i 
runs over all electrons in the 4/ shell. Parameters A nm 
do not depend on the ion, but are specific to the host 
compound (local environment of the ion) . Together with 
the ion-dependent parameters p n they describe the effect 
of the host crystal lattice on the rare-earth ion. Authors 
of paper j23| performed consistent analysis of the optical 
spectra of nine different rare-earth ions when they dope 
lanthanum trifluoride by substituting lanthanum ions, 
and deduced the values of parameters A nm for this com- 
pound. In our analysis, we neglect the possible difference 
between the local field in lanthanum trifluoride doped by 
rare-earth ions and the crystal field in generic rare-earth 
trifluoride compounds, and use the crystal field parame- 
ters derived in paper [231 ]. 

According to Eq. (jTU|l . in order to calculate (g\HpNc\g) 
we need to find the expectation value of [L x S]. To 
be specific, let us consider non-Kramers ions, i.e., ions 
that contain even number of electrons. In this case the 
ground state is non-degenerate and, as the direct calcu- 
lation shows, the expectation value ([L x S]) vanishes in 
the zero magnetic field. So, (-ffpNc) = 0, if B = 0, and, 
certainly, this agrees with the general statement that the 
usual weak interaction cannot generate linear Stark effect 
0- In the nonzero magnetic field the expectation value 
is of the form 

( 9 |[LxS]|p) = a[Bxk](B-k) 

+ 6[kB 2 -B(B-k)] , (14) 

where k is the unit vector orthogonal to the plane of 
hexagonal symmetry, and a and 6 are some constants de- 
rived from calculation. The first term in (|14|) corresponds 
to the simple structure that follows from most general ar- 
guments, as discussed in the beginning of this Section; a 
similar structure was considered in the literature 201 for 



Cs trapped in solid 4 He. The nature of the second term 
is not as obvious. An interesting feature of this term is 
that the vector 6k is a pseudovector: while a true vector 
is forbidden in centrosymmetric lattice, a pseudovector is 
generally allowed. We have confirmed numerically that 
the contribution proportional to 6k appears only due to 
4th and 6th multipolarities of the crystal field — indeed, 
one can construct a pseudovector from the irreducible 
tensors of 4th and 6th rank, but not from a 2nd rank 
tensor. For example, one can construct a pseudovector 
from the 4th rank irreducible tensor t, which represents 
the crystal field, in the following way: 

bki (X Cijk tjlmn ^Imnp ^rsuv ^uvpq ^qrsk • (1^) 

Here e^fc is the usual totally antisymmetric tensor, and 
summation over the repeated indices is implied. At zero 
temperature, the combination (|15l) does not appear ex- 
plicitly in the expression for the coefficient 6 — in this case 
the coefficient does not have a perturbation theory struc- 
ture, and, hence, the answer cannot be presented as a 
simple series in powers of crystal field. However, we also 
considered a nonzero temperature case. At high tem- 
perature, kT H c {,Als , the pseudovector 6k can be 
represented as a series in powers of H c {/kT. In this case 
numerics clearly indicate that, in agreement with (TSJ, 
the leading contribution to 6k appears in the fifth order 
in crystal field. 

The system of electron spins in RF3 at sufficiently low 
temperature undergoes transition to the ordered state. It 
is clear that the lowest temperature in the experiment is 
limited to the transition temperature. However, the ex- 
isting experimental data on RF3 is, unfortunately, scarce. 
It is known that the electro n sp ins in TbF3 order ferro- 
magnetically below 3.95 K j24j, and the spins in H0F3 
order antiferromagnetically at 0.53 K j^. From the re- 
sults of Ref. [2^| we also know that (Ce,Nd,Pr)F3 com- 
pounds remain paramagnetic at least down to 2 K. It is 
reasonable to conclude that other rare-earth trifluorides 
have transition temperature of ^ 1 K as well, so, in our 
estimates, we assume the temperature of 1 K. For non- 
Kramers ions where the ground state is not degenerate 
the effect is, actually, temperature independent as soon 
as the temperature in the compound is above the tem- 
perature of collective spin ordering ~ 1 K, and below the 
first energy interval in the crystal field splitting (~ 10 K 
for compounds where effect is the largest). However, the 
temperature dependence is very important for Kramers 
ions where the ground state is degenerate, and the effect 
there scales as 1/T. In the calculation we set T = 1 K. 

If we consider the external magnetic field of the order 
B < 1 T (tesla), the system is in the regime /j,bB <C 
kT <C Ae c f < Ae^g. Numerical analysis of the Hamil- 
tonian l|ll|) in this case gives the following parametric 
dependencies for the coefficients a and 6 describing the 
effect Q : 

non-Kramers ions: a oc — — - — , 6 oc 2 ; (16) 

BtfA L s A LS 
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0.36 


■10" 


-6 


0.27 
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0.93 
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-9 


0.13 


■10~ 5 


0.24 
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-5 
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TABLE II: Values of the coefficients a and b, Eq. (1141 . for different rare-earth ions. We assume that the magnetic field B in 
Eq. (I14H is expressed in units of T (tesla). The values of a and b for Kramers ions correspond to the temperature T = 1 K. 



Kramers ions: 



H, 



c£ 



kTA 



LS 



kT A 2 LS 



(17) 



The above expressions describe variation of the effect for 
all rare-earth ions, with the exception of Eu 3+ and Gd 3+ . 
Eu 3+ ion has the ground state with J = 0, so that the 
first energy interval is determined by the fine structure 
splitting rather than crystal field; in this case the effect 
includes additional orders of i? c fMLS : a °^ "ttS b oc 

A LS 

-jifS and so, is suppressed. The Gd 3+ ion has a half- filled 

4/-orbital with L = 0, so the expression ([L x S]), as we 
have already mentioned, vanishes completely for this ion. 
Table Inl lists the values of a and b for the rare-earth ions 
Ce through Yb derived from numerical calculations. 

The general features of the variation of a and b for dif- 
ferent ions can be understood with the help of Eqs. 
and 1)17(1 . In particular, the largest values of a among the 
non-Kramers ions correspond to the ions with the small- 
est first energy interval in the crystal field splitting. This 
agrees well with the formula 1(16(1 . Tb 3+ and Ho 3+ which 
have particularly large values of a, have very small first 
energy intervals ~ 5 cm -1 , whereas the average energy 
difference in the crystal field spectrum is ~ 100 cm"" 1 . 
Another general tendency observed for all ions is that 
the a-effect is larger than the 6-effect by approximately 
the factor H c {/Als, just as it should be according to 
Eqs. JH3), (JUJ. Although the expressions 1(16(1 and 1(17(1 
describe most general trends observed in Table ITT1 there 
is still considerable dispersion in the values. This is be- 
cause the many-electron states of the different ions have 
significantly different structure, and this kind of variation 
cannot be accounted for in a simple way. 

The rare earth ions in RF3 occupy the sites with local 
C2 symmetry. In addition, there is the Csj symmetry 
which describes different ionic sites in the lattice: for 
any arbitrary site, there are also sites rotated on 27r/3 
and —2tt/3 (rotation axis is the same as the local C2 
axis), as well as the sites with the inverted environment. 
The local site symmetry describes the geometry of the 
single site environment: it poses constraints on the crys- 
tal field parameter values, and hence, on any kinematics 
that appear locally. The existence of local ionic sites with 
different orientations in the lattice cell and the associated 
symmetry pose limitations on the form of any material 
(macroscopic) tensor. The kinematical relation 1(14(1 cor- 



responds to the macroscopically averaged quantity, and 
as such, is consistent with the combined Cqa symmetry 
(the local C%, plus C^a which describes the site orienta- 
tions). The values presented in Table[n] also, correspond 
to the kinematics averaged over the lattice. 

From equations ((10(1 and ((141) we derive the energy 
density induced by the weak interaction 

-Epnc = n^l-ffpNcIS') 

= anAE Qw(E-[Bxk])(B-k) , (18) 

the number density of rare-earth ions in the compound 
is about n ~ 1.8 • 10 22 cm" 3 . In the above equation 
we neglect the second term on the right-hand side part 
of Eq. lfH)l : although appearance of the pseudovector 
6k is itself an interesting point, in practice, as it follows 
from Table Ull the contribution of the 6-term is negligible 
compared to the contribution of the a-term. 

The macroscopic magnetization corresponding to 
Eq. CU) is 



M = - 



PNC 



dB 



iAE Q w |k(B • [k x E]) 

+ [kxE](B-k)} . (19) 



Numerical values of the coefficients a and A are presented 
in Tables [I] |nj we assume that the electric field is ex- 
pressed in units of V/cm, and the magnetic field B and 
magnetization M are expressed in units of T (tesla) . The 
temperature is set to T = 1 K, although, this is im- 
portant for Kramers ions only. From Tables [D] and 
the Eq. I(19|) we see that the effect is maximal for Dy 3+ , 
which is a Kramers ion, and Ho 3+ which is non-Kramers. 
Under favorable experimental conditions, it is reasonable 
to assume the applied magnetic field of B = 1 T and 
the electric field of E = 10 kV/cm. Then, the magnetic 
induction in the sample is 



DyF 3 : 4irM - 0.5 • 10~ n G 



HoF,: 4ttM ~ 0.1 • 10" 



G 



(20) 



Quadratic dependence of -Epnc on the value of mag- 
netic field seen in Eq. ((18|) is valid only if ^bB <§; kT for 
Kramers ions, and /is Ae c f for non-Kramers ions. For 
higher values of magnetic field the dependence turns to 
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linear. This means that for Dy and Ho compounds the 
PNC magnetization (|19fl increases linearly with B up to 
B~ 1 T and then stops to grow. 



IV. DISCUSSION 

The parity- violating weak interaction of electrons with 
the weak charge of rare-earth nuclei in rare-earth triflu- 
orides leads to the parity nonconserving magnetization 
presented in Eq. (|19|l . The effect depends on the exter- 
nal electric and magnetic fields. According to our esti- 
mates, the effect is maximal in Ho and Dy trifluorides. 
Assuming that the external magnetic field is B — 1 T, 
the external electric field is E = 10 kV/cm, and the tem- 
perature is T = 1 K, we find that the PNC magnetization 
is about 1CP 11 G. We can compare this estimate with the 
current experimental sensitivity that is under discussion 
in the literature. According to Ref. 0], it is possible to 
achieve the sensitivity of 3 • 10~ 16 G for 10 days of aver- 
aging with SQUID magnetometry, and it even might be 
possible to do 2 orders of magnitude better with appli- 
cation of different techniques. This leads to conclusion 
that the value of the PNC magnetization can, at least in 
principle, be measured to very good precision. 

However, an accurate theoretical treatment of the 
problem with precision of even several percents does not 
seem possible. At best, we can expect the accuracy at the 
level of 20-50%. Thus, extraction of the precise value of 
the weak charge directly from the results of experimen- 
tal measurements does not seem a viable possibility. It is 
possible, however, to look for the variation of the effect 
in a chain of nuclear isotopes in the same compound. 
In this respect, dysprosium, which has seven naturally 
occurring isotopes, looks most promising. This opens a 
possibility for the precise measurements of the Weinberg 
angle. There are several issues related to the hypotheti- 
cal measurements. The first is the normalization of the 
effect. Imagine, there are two samples made of different 
isotopes. They will be of slightly different shapes and vol- 
ume, and then one has to perform precise measurements 
of the usual electric and magnetic susceptibilities to nor- 
malize the PNC effect per unit volume. Another issue is 
related to the possible systematics. At the moment, we 
do not see any macroscopic systematic effects that would 
imitate the PNC magnetization ljlTi|) . Although Eq. (jT§|l , 
in essence, describes the nonlinear magnetoelectric effect, 
the crucial feature is that it violates parity at fundamen- 
tal level. In the lattice with center of inversion, the usual 
magnetoelectric effects cannot generate this correlation. 
A possible source of systematics can be due to the small 
local imperfections in the lattice which destroy the inver- 
sion symmetry. However such effects are of local nature 
and should average out to zero in the macroscopic sample 
with about 10 23 rare-earth sites. 

Another serious issue (H^E^) is the unknown value 
of nuclear neutron form factor which appears in the for- 
mula for effective nuclear weak charge — the quantity that 



is measured in atomic experiments. Electronic matrix 
element of the operator of weak charge interaction 
contains the expression p(r)Qw — p(r)(—N + Z(l — 
4 sin 2 ^)) averaged over the electron radial wave func- 
tions, p(r) here is the nuclear density which replaces the 
Dirac delta function in Eq. . It is known, however, 
that the density distribution is different for neutrons 
and protons, so the correct expression to be averaged 
is (—Np n (r) + Zp p (r)(l — 4sin 2 9w))i where p n {r) and 
p p (r) are the neutron and the proton nuclear density re- 
spectively. Thus, in a real experiment one probes the 
effective weak charge 

Q w = -Nq n + Zq p (l-4sm 2 e w ) (21) 

which contains the neutron and proton form factors q n 
and q p . The form factor q p can, at least in principle, 
be measured reliably with electron scattering, in muonic 
atoms, etc. It is expressed in terms of the proton distri- 
bution root-mean-square (RMS) radius R p , which is also 
the RMS radius of nuclear electric charge distribution. 
We assume that the uncertainty in the value of R p is 
negligible — the neutron RMS radius R n and the neutron 
form factor q n are known with much lower precision. 

Let us estimate the uncertainty in sin 2 9w due to un- 
certainty in the neutron form factor q n . To be specific, 
let us consider the suggested experiment with dyspro- 
sium trifluoride. There are seven stable isotopes of Dy: 
156, 158, 160, 161, 162, 163, 164. The odd and the even 
isotopes should not be considered together, because the 
odd isotopes then would give rise to nonmonotonic de- 
pendences due to the presence of the unpaired neutron. 
Clearly, it is best to analyze the even isotopes only. This 
leaves dysprosium-156, 158, 160, 162 and 164 — these are 
completely paired nuclei with practically the same de- 
formations, so the dependence of q n on the number of 
neutrons must be smooth and monotonic. 

Strictly speaking, in order to find the neutron form 
factor q n and the corresponding uncertainty one has to 
perform calculations of q n for deformed nucleus in the 
rotational s-wave state. However, to estimate the un- 
certainty, it is sufficient to use the known formula for 
spherical nucleus in the "sharp-edge" approximation |2?| 



1 



(22) 



The difference AR np = R n — R p was found in Ref. [23 
from analysis of experiments with antiprotonic atoms. 
The empirical fit in Ref. |2£j reads: 



AR np [fm] = (-0.04±0.03) + (1.01±0.15) 



N - Z 



(23) 



N + Z 

According to this expression the uncertainty in R n is 
5R n ~ 0.058 fm. From Ea. 122(1 . uncertainty in the neu- 
tron form factor is related to SR n in the following way 



Sq„ 



3 -(aZf^5R n 



(24) 
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The experiment with the chain of isotopes in the same 
compound would be measuring the ratio of Qw <|21|) for 
two different isotopes with A = Z + N and A' = Z + N': 

Qw -Ng n + Zq p (l - 4s 2 ) 

Q' w + Zq>(l - 4s 2 ) ' 1 > 

and thus, it would allow to eliminate the solid-state pre- 
factor and the associated theoretical uncertainty. Assum- 
ing there is no contribution from experimental error, the 
uncertainty in s 2 = sin 2 9w that corresponds to the error 
bounds in the Eq. (|23|l is: 

5s 2 ~ 5 • 1(T 4 . (26) 

Note, that the uncertainties in q n and q' n in Eq. I|25|l 
are correlated according to Eqs. lO- Moreover, 

main contribution to Ss 2 comes from the second term 
in Eq. I|23|) . while contribution of the constant term is 
negligible. 

There is also the experimental uncertainty SQw/Qw- 
In the possible experiment that measures the ratio 
Qw/Qwi this uncertainty translates into the uncertainty 



8s 2 in the value of sin 2 9w according to the formula: 

4ZAN [ Qw ) [ Qw ) 

This corresponds to the measurement with two dyspro- 
sium isotopes, 156 Dy and 164 Dy, AiV = 8. In order to 
match the uncertainty l]2(:iH . the experiment should be 
carried out with the precision exceeding SQw/Qw ~ 
8 • 1CP 5 . According to the estimates, statistical consid- 
erations allow for precision of the order ~ 3 • 10~ 5 , or 
better, depending on the magnetometry technique. 

The current experimental precision for the value of 
sin 2 6w measured at Z-peak, according to Ref. pj , is 
Ss 2 ~ 1.5 • 10 -4 . The solid-state experiment proposed 
in the present paper could achieve comparable level of 
precision, and, what is more important, would measure 
the value of the Weinberg angle in the domain of low 
energies (effective momentum transfer ~ 30 MeV). 

The authors would like to thank J. M. Cadogan and 
A. I. Milstein for stimulating discussions, and D. Budker 
and V. A. Dzuba for important comments. 
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